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Abstract 

A fully general relativistic description of the pulsar magnetosphere is provided. To be more 
concrete, a study of the pulsar magnetosphere is performed in the context of general relativistic 
magnetohydrodynamics (MHD) employing the so-called Grad-Shafranov approach. Not surpris- 
ingly, the resulting Grad-Shafranov equations and all the other related general relativistic MHD 
equations turn out to take essentially the same structures as those for the (rotating) black hole 
magnetosphere. Other different natures between the two cases including the structure of singular 
surfaces of MHD flows in each magnetosphere are essentially encoded in the different spacetime 
(metric) contents. In this way, the pulsar and the black hole magnetospheres can be described in 
an unified fashion. Particularly, the direction of poloidal currents circulating in the neutron star 
magnetosphere turns out to be the same as that of currents circulating in the black hole mag- 
netosphere which, in turn, leads to the pulsar and the black hole spin-downs via the "magnetic 
braking" . 
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I. Introduction 



Radio pulsars are perhaps tlie oldest-known and the least energetic among all of the 
pulsar categories. The thoretical study of these radio pulsars can be traced back to the 1969 
work of Goldreich and Julian In their pioneering work, Goldreich and Julian argued 
that pulsars, which are thought to be the rotating magnetized neutron stars, must have a 
magnetosphere with charge-separated plasma. They then demonstrated that an electric force 
which is much stronger than the gravitational force will be set up along the magnetic field 
and as a result, the surface charge layer cannot be in dynamical equilibrium. Then there 
appear steady current flows along the magnetic field lines which are taken to be uniformly 
rotating since they are firmly rooted in the crystalline crust of the pulsar surface. Although 
it was rather implicit in their original work, this model for the pulsar electrodynamics does 
indeed suggest that the luminosity of radio pulsars is due to the loss of their rotational 
energy (namely the spin-down) and its subsequent conversion into charged particle emission 
which eventually generates the radiation in the far zone. 

Since the pioneering proposals by Gold [3| and by Pacini P|, it has by now been widely 
accepted that indeed pulsars might be rotating magnetized neutron stars. Nevertheless, since 
then nearly all the studies on the pulsar electrodynamics have been performed by simply 
treating the region surrounding the rotating neutron stars as being fiat. This simplification 
may be sufficient just to gain some insight into rough understanding of the origin of the 
pulsars' radiation. However, nearly over thirty years have passed of the study of pulsar 
electrodynamics and it seems to be that now we should treat the problem in a more careful 
and rigorous manner, namely, in a fully relativistic fashion. Here by relativistic treatment, 
we mean that we are dealing with the highly self-gravitating compact objects and of course 
we shall have the rotating neutron stars in mind. Thus we first begin by providing the 
rationale for treating the vicinity of spinning compact objects such as a rotating neutron star 
relativistically as a non-trivial curved spacetime. Indeed^e now have a considerable amount 

h n 

of observed data for various species from radio pulsars p| to (anomalous) X-ray pulsars |^ . 
Even if we take the oldest-known radio pulsars for example, it is not hard to realize that these 
objects are compact enough which really should be treated in general relativistic manner. To 
be more concrete, note that the values of the parameters characterizing typical radio pulsars 
are ; ro(radius) ~ 10®(cm), M(mass) ~ IAMq ~ 2 x 10^^{g), r(pulsation period) ~ 10"'^ — 
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l(sec), _B(magnetic field strength) ~ lO^^(G'). Thus the Schwarzschild radius (gravitational 
radius) of a typical radio pulsar is estimated to be rsch = 2GMq/c^ ~ 3 x lO^(cm) (where 
G and c denote the Newton's constant and the speed of light respectively) and hence one 
ends up with the ratio VQ/rsch ~ 10®(cm)/3 x lO^(cm) ~ 3. This simple argument indicates 
that indeed even the radio pulsar (which is perhaps the least energetic among all of its 
species) is a highly self-gravitating compact object that needs to be treated relativistically. 
Once again, the "relativistic" treatment here means that the region surrounding the pulsar, 
namely the pulsar magnetosphere has to be described by a curved spacetime rather than 
simply a flat one. Then the question now boils down to ; what would be the relevant 
metric to describe the vicinity of a rotating neutron star ? Although it does not seem to 
be well-known, fortunately we have such a metric of the region exterior to slowly-rotating 
relativistic stars such as neutron stars, white dwarfs and supermassive stars and it is the 
one constructed long ago by Hartle and Thorne j|5|. Thus the Hartle-Thorne metric is a 
stationary axisymmetric solution to the vacuum Einstein equation and in the present work, 
we shall take this Hartle-Thorne metric as a relevant one to represent the spacetime exterior 
to slowly-rotating neutron stars. 

Having provided the rationale for treating the region surrounding a magnetized rotating 
neutron star, namely a pulsar fully relativistically as a curved spacetime, we now state 
our particular objective in this work. In the present work, we would like to study the 
pulsar magnetosphere in the context of general relativistic magnetohydrodynamics (MHD) 
by employing the so-called Grad-Shafranov approach We shall consider both the force- 
free and full MHD situations and accordingly derive the Grad-Shafranov equations for each 
case, namely, the pulsar equation and the pulsar jet equation respectively. Not surprisingly, 
then, the resulting Grad-Shafranov equations and all the other related force-free equations 
or general relativistic MHD equations turn out to take essentially the same structures as 
those for the (rotating) black hole magnetosphere . The essential distinction between 

the two cases, however, is the spacetime (metric) contents. For the pulsar magnetosphere 
case, one needs to choose the Hartle-Thorne metric mentioned above whereas for the black 
hole „„phe.e ease, one to select the Ke„ blaek hole .netnc Q. Then as a 
consequence of this, the singular surfaces of MHD flows in the pulsar magnetosphere and 
those in the black hole magnetosphere exhibit substantially different nature. In this way, the 
pulsar and the black hole magnetospheres can be described in an unified fashion. This last 
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point, namely, an unified picture of both tlie pulsar and the black hole electrodynamics is the 
key proposal of the present work. There is, however, as strong motive for the present work 
and it is the uncomfortable current state of affair that there still is no generally accepted 
standpoint about the structure of pulsar magnetosphere yet To be more concrete, there 
is yet no complete model for the structure of longitudinal (or poloidal) currents circulating 
in the neutron star magnetosphere that can provide the solution to the problem, say, of 
pulsar spin-down (namely, the "braking" of rotating neutron star). As we shall see shortly 
in the text, a partly satisfying solution to this problem will be provided by treating the 
region outside a magnetized rotating neutron star as a curved spacetime represented by the 
Hartle-Thorne metric. Namely it turns out that the structure of space charge-separation 
and the direction of poloidal current particularly in the force-free limit correctly lead to 
the magnetic braking torque that spins down the rotating neutron star. This implies that, 
particularly in the technical aspect, the mechanism of Blandford-Znajek type can be adopted 
to provide a model which takes the magnetized rotating neutron star as the central engine 
for some radio. X-ray and even (soft) gamma-ray astrophysical phenomena ^ just as it has 
been employed to construct a model taking the rotating black hole as the central engine 
for active galactic nuclei (AGNs)/quasars 12] or even gamma-ray bursts (GRBs) 13]. It, 
however, seems fair to say that this observation should not be taken as a new discovery but as 
something that has been expected to some extent. And it is due to the fact that historically 
the Blandford-Znajek mechanism has been strongly motivated by and thus constructed from 
the original pulsar model of Goldreich and Julian [l| (and perhaps others) with the purpose 
to apply the main idea to the case of rotating black hole. We also note that there have 
been extensive critical examinations of the operational aspect of the Blandford-Znajek type 



mechanism carried out by Punsly 



II. Electrodynamics around the slowly-rotating neutron stars 

1. The Hartle-Thorne metric for the region exterior to the slowly-rotating 
neutron stars 

The Hartle-Thorne metric ^ is given, in terms of the (3 + l)-split form, by 

ds'^ = -a^dt^ + grrdr"^ + gBede"^ + 5-00(^0 + 13'^dtf (1) 
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where the lapse a, (angular) shift (3'^ and the metric components are 

2J 



a 

9rr 

9tt 



AR, p^ = -uj = — 



S_ 
A' 



VD^ = rMsin^6', 



(2) 



■[AR 



Asin^ 6], 



r 



and 



A 
R 

S 
A 



2M 2J2 



1 + 2 
1 - 2 
1 + 2 



r J' 



1 + 



+ 



M\ 5Q-J^/M 

8 M3 
5M\ 5Q-J^/M 



Ql 



+ ^ 



. r 

r 



1) P2(C0S( 



(3) 



1 + 



r 

2M\ 



8 M3 



Q2(^-l) ^2(cos^) 



5 Q - JVM 
^8 M3 



2M 



rQ^(^-l)-Q2(^-l)} A(cos( 



[r2(l-2M/r)]V2 

with M, J and Q being the mass, the angular momentum and the mass quadrupole moment 
of the (slowly) rotating neutron star respectively, P2(cos^) = (3cos^^ — l)/2 being the 
Legendre polynomial, and being the associated Legendre polynomial, namely, 
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As is well-known, the only known exact metric solution exterior to a rotating object is the 
Kerr metric IKI] . Thus it would be worth clarifying the relation of the Hartle-Thorne metric 
or slowly-rotating relativistic stars given above to the Kerr metric. As Hartle and Thorne 
^ pointed out, take the Kerr metric given in Boyer-Lindquist coordinate and expand it 
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to second order in angular velocity (namely, the angular shift f3'^) followed by a coordinate 
transformation in the (r, 6')-sector, 



^ - cos ^ sin 0^7(1 + —) (6) 

where a = J/M. Then one can realize that the resulting expanded Kerr metric coincides 
with the particular case Q = J^/M (with Q being the mass quadrupole moment of the ro- 
tating object) of the Hartle-Thorne metric. Therefore, in general this Hartle-Thorne metric 
is not a slow-rotation limit of Kerr metric. Rather, the slow-rotation limit of Kerr metric 
is a special case of this more general Hartle-Thorne metric. As a result, the Hartle-Thorne 
metric with an arbitrary value of the mass quadrupole moment Q can generally describe a 
(slowly- rotating) neutron star of any shape (as long as it retains the axisymmetry) . 
Next, since we shall employ in the present work the Hartle-Thorne metric to represent the 
spacetime exterior to slowly-rotating neutron stars, we would like to carefully distinguish 
between the horizon radius of the Hartle-Thorne metric and the actual radius of the neutron 
stars. Indeed, one of the obvious differences between the black hole case and the neutron star 
case is the fact that the black hole is characterized by its event horizon while the neutron 
star has a hard surface. Since this solid surface of a neutron star (which we shall henceforth 
denote by tq) lies outside of its gravitational radius which amounts to the Killing horizon 
radius of the Hartle-Thorne metric, th at which A = (1 — 2M/r + 2J^/r^) = in eq.(3), we 
have Tq > r^- In the present work, however, we shall never speak of the Killing horizon of 
the Hartle-Thorne metric as it is an irrelevant quantity playing no physical role. 
Now we turn to the choice of an orthonormal tetrad frame. And we shall particularly 

n 

choose the Zero-Angular-Momentum-Observer (ZAMO) \l^ frame which is a sort of fidu- 
cial observer (FIDO) frame. Generally speaking, in order to represent a given background 
geometry, one needs to first choose a coordinate system in which the metric is to be given 
and next, in order to obtain physical components of a tensor (such as the electric and mag- 
netic field values), one has to select a tetrad frame (in a given coordinate system) to which 
the tensor components are to be projected. As is well-known, the orthonormal tetrad is 
a set of four mutually orthogonal unit vectors at each point in a given spacetime which 
give the directions of the four axes of locally-Minkowskian coordinate system. Such an or- 
thonormal tetrad associated with the Hartle-Thorne metric given above may be chosen as 



adt = {/^Rf/'^dt, 
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The local, stationary observer at rest in this orthonormal tetrad frame e"^ has the worldline 
given by {dr = 0, d9 = 0, (dcj) + [3^dt) = 0} which is orthogonal to spacelike hypersurfaces 
and has orbital angular velocity given by 

n. , 1 

(9) 



dt 
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This is the long-known Lense-Thirring precession [ly] angular velocity arising due to the 
"dragging of inertial frame" effect of a stationary axisymmetric spacetime. Indeed, it is 
straightforward to demonstrate that this orthonormal tetrad observer can be identified with 
a ZAMO carrying zero intrinsic angular momentum with it. To this end, recall that when 
a spacetime metric possesses a rotational (azimuthal) isometry, there exists associated ro- 
tational Killing field = [d/dcf))^ = 6^ such that the inner product of it with the tangent 
(velocity) vector u'^ = dx^/dr (with r denoting the particle's proper time) of the geodesic 
of a test particle is constant along the geodesic, i.e.. 



gapnfu'^ = g^tm'*'u^ + g^^m'*'u^ 
dt d(f) 



(10) 



Now, particularly when the local, stationary observer, which here is taken to be a test 
particle, carries zero angular momentum, L = 0, its angular velocity becomes 

d(j) {d(j)/ dr) gt4> 



dt (dt/dr) 



(11) 



and this confirms the identification of the local observer at rest in this orthonormal tetrad 
frame given above with a ZAMO. 

2. Electrodynamics in curved spacetime : the (3 + l)-split formalism 

Generally speaking, when dealing with the electrodynamics in a curved spacetime, rela- 
tivists prefer a geometric, covariant, frame-independent approach, representing, say, the 
electromagnetic field by the field strength tensor F^y. The astrophysicists, on the other 
hand, would prefer to split this tensor into a 3-dimensional electric field E and magnetic 
field B, sacrificing the general covariance of the theory in order to get some insight and 
achieve a comparison with the familiar flat spacetime electrodynamics. As has been first 
developed by Macdonald and Thorne ^ , fortunately in stationary curved spacetimes, such 
as those outside the rotating black holes (the Kerr metric) and the rotating neutron stars 
(the Hartle-Thorne metric discussed above), one can actually reformulate electrodynamics 
in terms of an absolute 3-dimensional space and an universal time. And the variables in 
this reformulation are the familiar electric and magnetic fields (E,B), charge and current 
density (pe,j)- Indeed, this absolute-space/universal-time formulation of stationary curved 
spacetimes has deep roots in the so-called (3 + l)-split formalism of general relativity which 
had originally been employed in the canonical (or Hamiltonian) quantization of gravity in 
the 1950s and is nowadays being used in the numerical relativity. Normally, however, such 
(3 + l)-split formalism has not been welcome by most relativists due to the arbitrariness 
of the choice of fiducial reference frame. In the case of stationary black hole/neutron star 
electrodynamics, however, there is one set of fiducial observers preferred over all others : 
the Zero-Angular-Momentum-Observer (ZAMO) or Locally-Non-Rotating-Frame (LNRF) 
observer that we discussed in the above subsection. As is well-known and as we shall see 
in a moment, when one uses this ZAMO frames one realizes that the (3 + 1) equations of 
black hole/neutron star electrodynamics are nearly identical to their counterparts in fiat 
spacetime electrodynamics. For general formulations and more detailed discussions of this 
(3 + l)-split formalism we refer the reader to jj]. As such, throughout this work, we shall 
also employ this space-plus-time formalism in which all the physical quantities are repre- 
sented by 3-dimensional scalars and vectors as measured by ZAMO, the local observer. For 
instance, the physical electric and magnetic field components as measured by ZAMO can 
be read off as the projection of F^j^y onto the ZAMO orthonormal tetrad frame (eq.(8)). 
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Fab = F^y{e^AeB) and Fab = {Fio, Fij} where 

Ei = FiQ, E = {Ei}, 

= 2^ijkF^^, B = {Bi}. 

III. Pulsar equation - The force-free limit of the Grad-Shafranov equation 

The so-called Pulsar equation refers to the force-free limit of the more general Grad- 
Shafranov equation. The force-free condition essentially amounts to ignoring the (inertial) 
contributions of the plasma particles when the plasma energy density is assumed to be sub- 
stantially smaller than that of the magnetic field. 
1. Basic equations 
(1) Force- free condition 

In order eventually to describe the force-free pulsar magnetosphere, we start with the 
Maxwell equations in the background of the stationary axisymmetric rotating neutron star 
spacetime |7| 

V ■ E = 47rpe, V ■ B = 0, 

V X (aE) = (B ■ Vcj)m, (12) 

V X (aB) = 47raj - (E ■ Vtu)m 

where we dropped the terms £k(---) = 0, £m(---) = due to stationarity and axisymmetry. 
And here k = [d/dt) and m = {d/dcj)) denote the time-translational and the rotational 
Killing fields associated with the stationarity and the axisymmetry of the Hartle-Thorne 
metric, respectively and hence m ■ m = g^^ = and m = roe^ = {g(j,(f,Y^'^e^. Note also 
that since all the measurements are made by ZAMO, the lapse function a is introduced to 
convert the ZAMO's proper time dr over to the global time dt. Throughout in this section, 
the force-free condition is assumed to hold, i.e., 

PeE+jxB = 0, B = Bt + Bp (13) 

which also implies the degeneracy condition E ■ B = 0. Then this force-free condition 
indicates that the charged particles are fiowing along the magnetic field lines and hence the 
toroidal (angular) velocity of magnetic field lines (which are frozen into plasma) relative to 
ZAMO is given by 
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and then j — jr + jp with = Pe^r where Vt consists of Vp given above and the streaming 
velocity along the toroidal magnetic field fines. Then from the force-free condition above, it 
follows that 

E = Ep = -— Jt X Bp. (15) 

Pe 

(2) Poloidal field components 

Consider a magnetic flux through an area A whose boundary is a m-loop, 
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(3) Poloidal current and toroidal magnetic field 

Consider a poloidal current through the same area whose boundary is a m-loop, 
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then similarly to the case of poloidal magnetic fleld, V/ = —2Trzue^ x (ctjp), we get 



and then from eqs(17), (20), 
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Next, consider the Ampere's law in Maxwell equations 
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Next, using eqs.(12) and (18), one gets the Gauss law equation 
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while the Ampere's law in eq.(12) yields 
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Alternatively, by combining eqs.(24) and (25), one gets 
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2. The Grad-Shafranov approach 

(1) The Grad-Shafranov equation 

Consider the force-free condition peE + j x B = 0, and focus on its "poloidal component" 
equation, 

-pe(vF X B)|p+j X B|p = 

which yields. 
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Now by plugging eqs.(24) and (25) in (27) above and using eq.(23), one arrives at 
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where we assumed a general situation in which the magnetic field lines are, although rooted 
in the pulsar surface, allowed to have differential rotation. This is the stream equation to 
determine the field structure of the force-free pulsar magnetosphere. 
(2) Electric charge and toroidal current density 
rom eqs.(26) and (27), one gets the expressions for the charge and toroidal current density 
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(3) Singular surfaces 

In order to distinguish the singular surfaces in the pulsar magnetosphere from those in the 

(rotating) black hole magnetosphere, we first note the generic distinction between the black 
hole magnetosphere and the pulsar magnetosphere. 
{Black hole magnetosphere) 

In this case, the field angular velocity is generally in no way connected with the angular 
velocity of the black hole culrn) — ^h- Indeed, {ftp — uj) changes sign from minus to plus 
as one moves away from the symmetry axis (recall that cu denotes the angular velocity of 
ZAMO). 

{Pulsar magnetosphere) 

In this case, all the magnetic field lines are firmely rooted in the crystalline crust of the 
pulsar surface, namely Jlp = Jljvs > cu. Thus {ftp — cu) > namely, since ftp — ^ns, the 
field angular velocity should be greater than that of ZAMO, lo, everywhere. 
To summarize, the spinning black hole rotates the "surrounding" magnetic field lines (gener- 
ated, say, by currents in the accretions disc) essentially via the frame dragging effect whereas 
the neutron star rotates "its own dipole" magnetic field lines via its spin motion itself. 
We start with the light cylinder. When treating the region exterior to the rotating neutron 
star as a fiat spacetime, there was a single light cylinder at uj — c/Qp- For the case at 
hand where the region outside of the rotating neutron star is described by a generic curved 
spacetime, there still is a single light cylinder where the denominators of pe and jr vanish, 
namely at 

etc 

ruL = r. 30 

[ilf — OJ) 

The only change from the fiat spacetime treatment to the curved spacetime one is the notion 
of relative angular velocity with respect to ZAMO since now all the measurements are made 
by a local fiducial observer which is ZAMO. Particularly note that in the present case of 
rotating neutron star, there is only one zero for the denominators of Pe and in eq.(29) 
instead of two since the field angular velocity should be greater than that of ZAMO, oo 
everywhere as we explained above. Certainly, this is in contrast to what happens in the 
case of rotating black hole magnetosphere when there are two light cylinders where the 
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denominators of pe and jx vanish, namely at 
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On the hght cyhnder, the numerators should vanish in order to have finite pe and jt there, 
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which is called the "critical condition" . 

We now point out the implication of this distinction between the single light cylinder in the 
pulsar magnetosphere and the double light cylinders in the black hole magnetosphere. 
In the case of rotating black hole magnetosphere, the source or origin of the plasma that 
will fill the magnetosphere is understood as follows. Note that as one moves from the inner 
light cylinder toward the outer one, the angular velocity of the magnetic field lines grows, 
namely {Qp < u) {Qp > uj). Thus somewhere between the two light cylinders, there 
should be the null surface where 
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and hence Eip — —Vp x Bp — 0. Then on this null surface. 
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and this null surface occurs ai w = wn {"^il < < ^ol)- Indeed we refer to it as 
the "null surface" since there must be the spark gaps or the creation zone of nearly neutral 
plasma situated in its neighborhood. The velocity of the magnetic field lines relative to 
ZAMO 

'Q.p — uj' 



\Yp\ = 
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begins to increase from zero at ci7 = cuat toward c at = vjol and then further to infinity 
as one moves far away from the outer light cylinder where a; — > 0, a ^ 1. As one moves 
inward, on the other hand, it begins to increase in magnitude (upon changing sign) from 
zero w — W]^ toward —c at zu — zuil and then again further to (negative) infinity at the 
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horizon r/^ as a = a{rH) = there. Thus ZAMOs in the outer region of the magnetosphere 
zo > wn see the centrifugal "magnetic shngshot wind" blowing outward from the vicinity of 
the null surface to the acceleration zone and ZAMOs in the inner region of the magnetosphere 
zu < zon see the centrifugal magnetic slingshot wind blowing inward to the horizon. 
To summarize, the global structure of the pulsar magnetosphere associated with the singular 
surfaces is indeed quite different from that of rotating black hole magnetosphere. And it 
indeed is related to the fact that the black hole is characterized by its event horizon while 
the neutron star has a hard surface. 

3. Poloidal current in the neutron star magnetosphere and the pulsar spin- 
down 

As advertized earlier in the introduction, we now address the issue of pulsar spin-down 
essentially due to the magnetic braking torque in terms of the structure of longitudinal (or 
poloidal) currents circulating in the neutron star magnetosphere. To this end, we should 
start with the space charge-separation in the force-free limit. The presentation that will 
be given below for the present case of pulsar (which is being treated in a fuUv general 
relativistic manner) essentially follows that for the case of rotating black hole jq| in the 
context of Blandford-Znajek mechanism. 

Recall the definitions for the poloidal magnetic flux (or stream function) \E' and for the 
poloidal current I given, respectively, by 

^ = [ B-dS, I = - [ aj-dS. 

J A J A 

First, we consider the case when the angular momentum J^s and the (asymptotic) direction 
of the magnetic field B are parallel. We begin by noting that the magnetic flux (and B) 
is defined to be positive when it directs upward while the poloidal current is defined to be 
positive when it directs downward. Normally one imposes the condition of no net loss of 
charge from the pulsar. This amounts to demanding that the net current flowing into and 
out of the neutron star surface A^s vanishes, namely 

/ «j ■ rfS = / (-§-] Bp.dS = - I{^o)] = (34) 

JAns JAns \ J 

where we used j = jr + jp, jr ■ c^S = and d"^ = B • = Bp ■ dS. And for the expression for 
the poloidal current density jp, we used eq.(21). Here, the surface integral is taken only over 
the (nothern) hemisphere of the neutron star surface from the (north) pole where \E' = \l'o 
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to the equator where \I/ = "^eq as the pulsar's intrinsic dipole moment would generate dipole 
magnetic fields. This condition indeed implies the presence of some "critical" magnetic 
surface = \l/c such that 

> (for ^0 < ^ < ^c) : Region I, 
^ = (for VP = vl/J, (35) 
< (for < ^ < ^eq) ■ Region II. 

Then one can realize from eqs.(35) and (21) that the electric current fiows inwardly along 
the magnetic field lines for \l/o < ^ < from the acceleration region to the neutron star 
surface while it fiows outwardly for < ^ < from the surface to the acceleration region. 
And on the critical magnetic surface \1/ = \l/c, the net current is zero (i.e., infiowing current 
= outfiowing current) but the charge density there may not be exactly zero and hence the 
critical magnetic surface may not serve as exactly a "charge-separating" surface. This point 
can be envisaged from the expression for the charge density pe in eq.(29) where one can 
notice that dl /d^ is zero but G ■ V\l/ is not necessarily zero along = Indeed, this last 



quantity can roughly be identified with t 
Bz). In the Goldreich- Julian (GJ) model 



the vertical component of the magnetic field (i.e., 
l| of purely charge-separated magnetosphere, the 
critical magnetic surface has been defined as the one on which i?^ = all the way. In our 
model of pulsar magnetosphere with infinite supply of quasi-neutral plasma, on the other 
hand, it is being defined as the one on which dl /d^ = 0. As such, generally the critical 
magnetic surfaces in the two models may not coincide and as a result, in our model the 
critical magnetic surface may not act as precisely a charge-separating surface. The origin of 
different predictions of these two models shall be discussed in detail later on in the subsection 
5 of the present section. 

Next, in the case when the angular momentum J^s and the (asymptotic) direction of the 
magnetic field B are antiparallel, all the quantities involved would carry fiipped signs, namely 

< (for ^0 < ^ < ^c) : Region I, 
^ = (for VP = vl/J, (36) 
> (for < ^ < ^eq) ■ Region II. 

This second case when J7V5 (or Jbh that we shall also discuss shortly) and B are antiparallel 
is indeed hkely to happen. First for the pulsar case, its spin angular momentum Jat^ and its 
intrinsic magnetic dipole moment may well be aligned and pointing the opposite directions 
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at the same time to yield the antiparallel configuration of this type. For the rotating black 
hole case, on the other hand, consider, for instance, that the poloidal magnetic fields come 
from the toroidal currents in the accretion disc around the hole. Clearly, the hole and the 
disc would be corotating but if the excess charge is due to that of ions, the spin of the hole 
Jbh and the magnetic field would be parallel whereas if it is due to that of electrons, the 
two would be antiparallel instead. 

With this preparation, we now turn to the determination of the structure of space charge- 
separation and the direction of the poloidal current in the neutron star (and the rotating 
black hole for comparison) magnetosphere. Consider the charge density in the neutron star 
magnetosphere given earlier in eq.(29) (note that its structure remains the same even in the 
rotating black hole magnetosphere except that the associated spacetime metric content is 
distinct between the two cases), 



Pe 



m 


f -G-Vvl/ 


a 







(37) 



Now using this expression for the charge density along with eqs.(35) and (36), one can, in 
principle determine the sign of the separated charges in different domains of the neutron star 
magnetosphere. The directions of poloidal current densities are determined using basically 
the eq.(21) accordingly. The resulting charge-separation and the poloidal current direction 
are depicted in FIG.l. It is indeed quite instructive to compare the present case of pulsar 
magnetosphere structure with that of rotating black hole magnetosphere structure. The 



latter had been studied in detail in the literature llTl l and here in the present 

work, we have elaborated on it by further considering the case when the spin of the hole 
and the (asymptotic) direction of the magnetic field are antiparallel - see FIG.2. 
In practice, however, determining the sign of the separated charges in different domains of 
the neutron star/black hole magnetosphere using eqs.(35), (36) and (37) is by no means a 
straightforward job. Technically, the associated difficulty arises from the fact that we need 
to figure out which term, between KStt"^! /a^){dl /d'$)\ and |G ■ V\l'| in the numerator of 
eq.(37), is greater than the other in different domains of the magnetosphere. By contrast, the 
determination of the structure of charge separation (i.e., the sign of separated charges) in the 
GJ model (but only inside the light cylinder) that we shall discuss shortly and summarized in 
FIG. 3 using eq.(38) below is rather straightforward. Fortunately, however, there is a guiding 
principle that allows us to determine the sign of separated charges in different domains with 
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^IG. 1: Pulsar magnetosphere - a general relativistic treatment. The presence of accretion disc 
2nl | here is not mandatory but is assumed for parallel comparison with the case of black hole 
magnetosphere. The star's intrinsic dipole moments are not drawn in the figures, 
(a) When Jns and B are parallel, (b) when Jns and B are antiparallel. 

confidence. And that is just the insightful realization that the different domains of the 
magnetosphere around the compact object such as the rotating neutron star or black hole 
should rotate in the same direction as the compact object itself ( "corotation" ) if they rotate 
faster than ZAMO, the local inertial observer carryng out all the observations, whereas they 
should rotate in the opposite direction ( "counter- rotation" ) provided they rotate slower than 
ZAMO. In order to work out this guiding principle and determine the structure of charge 
separation as a result, we begin by defining the angular momentum of the magnetosphere. 
The pulsar magnetosphere consists of the poloidal magnetic field generated essentially by the 
intrinsic dipole moment of the neutron star and the poloidal electric field generated both by 
the toroidal current-poloidal magnetic field in the force-free limit (see eqs.(14) and (18)) and 
by the separated space charge. And particularly these two sources of the poloidal electric 
field are closely related since the toroidal current arises as a result of the angular motion 
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of the poloidal magnetic field lines (being dragged along by the spin of the neutron star) 
along which the plasma flows (in the force- free limit). Therefore if one can determine the 
direction of the poloidal electric fleld, the structure of charge separation can be determined 
accordingly. Besides, since the motion of the plasma eventually generates the poloidal 
electric fleld, the angular momentum of the magnetosphere can be deflned in terms of its 
poloidal electromagnetic fleld. In general, the angular momentum of an electromagnetic fleld 
is deflned by Jem = J d^x[r x (Ep x Bp)]. Thus we only need to determine the direction 
of this poloidal electric fleld Ep (for a given poloidal dipole magnetic fleld Bp) that leads 
to either corotation or counter-rotation of the different domains of the magnetosphere (as 
measured by ZAMO, the local inertial observer) with the rotating neutron star/black hole, 
i.e., Jem ~ ±JAf5,_BH- This is how the structure of charge separation in different domains of 
the pulsar magnetosphere summarized in FIG.l and that of the black hole magnetosphere 
summarized in FIG2. have been actually flxed. It is worthy of note that for the case of 
pulsar magnetosphere, the sign of separated charges in different domains such determined 
is the same as that in the original GJ model that we shall turn to in a moment except that 
now for our force- free treatment, it can be extrapolated outside the light cylinder. Next, for 
the case of black hole magnetosphere, the structure of charge separation such determined 

n 

turns out to be exactly the same again as the result derived by Okamoto [8] some time ago 
via different reasoning. Then to summarize, it is not surprising (since it has been expected 
to some extent) but still interesting to realize that the structure of magnetosphere of the 
pulsar and the rotating black hole are not quite the same let alone the different structure of 
singular surfaces that we stressed earlier. If we emphasize it once again, this difference can 
be attributed to the fact that all the magnetic fleld lines are flrmely rooted in the crystalline 
crust of the pulsar surface and hence flp = ^ns > ^ namely, the fleld angular velocity is 
greater than that of ZAMO, u, everywhere. Consequently, to the ZAMO of rotating neutron 
star, the whole pulsar magnetosphere looks corotating and thus the resulting direction of 
the poloidal electric fleld eventually determines the structure of space charge separation as 
given in FIG.l. In the black hole case, however, the fleld angular velocity fip is generally in 
no way connected with the angular velocity of the black hole uj{rH) = Qh- Indeed, [Qp — uj) 
changes sign from minus to plus as one moves away from the symmetry axis (recall that u 
denotes the angular velocity of ZAMO). Namely, inside the "null surface" that we discussed 
earlier in subsection 2, ZAMO rotates faster than the magnetic fleld lines (and hence than 
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FIG. 2: Rotating black hole magnetosphere. Although the critical magnetic surfaces have been 
sketched more or less as straight lines in these figures, they would in fact be curved in actual 
situation. 

(a) When Jbh and B are parallel, (b) when Jbh and B are antiparallel. 

the inner part of the magnetosphere) while outside of it, ZAMO rotates slower than the 
magnetic field lines (and thus the outer part of the magnetosphere). As a result, to the 
ZAMO of rotating black hole, the part of black hole magnetosphere outside the null surface 
looks corotating whereas its part inside the null surface appears to counter-rotating. Prom 
this one can realize the directions of the poloidal electric fields which, in turn, determines 
the structure of space charge separation as given in FIG. 2. And for both pulsar and rotating 
black hole cases, it is rather straightforward to see that the structure of charge-separation 
and the direction of longitudinal (or poloidal) current (denoted in the figures by /) circulating 
in the magnetospheres actually lead to the magnetic braking torques, namely the Lorentz 
torques N = !ai^s{h)^^ ^ ^ ^^^^ ^Pi^ down the rotating neutron star and the 

black hole as it is always directed opposite to the spins regardless of whether the spin and 
the (asymptotic) direction of the magnetic field are parallel or ntiparallel. Here, the surface 
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area A^s{h) over which the integral of the non- vanishing Lorentz torque density is to be 
taken might need some careful clarification. For the case of rotating neutron star, this area 
obviously should be its surface where the poloidal current crosses the magnetic field lines 
and hence generates non-vanishing braking torque. For the case of black hole, however, 
the nature of this area might seem quite ambiguous but our suggestion here is to invoke 
ihe notion of "stretched horizon" as an incarnation of the so-called "membrane paradigm" 
18| . Indeed, the philosophy that underlies the membrane paradigm is an attempt to have 
an intuitive picture of Blandford-Znajek mechanism by first assuming the appearance of 
stretched horizon (just outside the event horizon) and then introducing (fictitious) surface 
charge and current density on it. Then one of the most intriguing consequences of such 
assumption is that if we choose to do so, the (stretched) horizon behaves as if it is a conductor 
with finite resistivity. To be more specific, since there are now both current and resistivity 
on the horizon, one might naturally wonder what would happen to the Joule heat generated 
when those surface currents work against the resistance and how it would be related to 
the electromagnetic ene rgy going down the hole through the horizon. Indeed, Znajek and 
independently Damour [1^ provided a simple and natural answer to this question. Namely, 
they showed in a consistent and elegant manner that the total electromagnetic energy flux 
(i.e., the Poynting flux) into the rotating Kerr hole through the horizon is indeed precisely 
the same as the amount of Joule heat produced by the surface currents when they work 
against the surface resistivity of in. Therefore, motivated and encouraged by these ground 
works for the advent of the membrane paradigm, here we also assume that the surface 
area over which the integral of the non- vanishing Lorentz torque is to be taken is just this 
stretched horizon with surface current. Then the real poloidal current in the black hole 
magnetosphere and this surface current on the horizon together are supposed to complete 
the circuit. 

Now to summarize, this unified picture can be thought of as a satisfying solution to both the 
magnetized rotating neutron star interpretation of radio/X-ray pulsars and the rotating 
supermassive black hole interpretation of AGNs/quasars 12] or even GRBs [l^ . 
It is, however, the following point that is of great interest and has been the strong motive for 
the present study. And it is the difference in the nature of the origin/source of charges which 
get separated in the domains of the magnetosphere and of the resulting poloidal currents 
between the two pulsar models - ours and that of Goldreich- Julian's. It goes as follows. 
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First, the original GJ model can be thought of as the purely charge-separated solution in 
which only one charge species can be assigned at a given point in space. As a result, the 
poloidal current in this charge-separated solution is directly proportional to charge density. 
Now, the origin/source of charges in this GJ model is basically the surface of the pulsar 
itself and the poloidal current flows only until these charges get separated (by the strong 
local electric field near the surface of the star) and then reach the equilibrium GJ charge 
density given by 

GJ ^ 5-17 r2 Bz /oo\ 

^"2^1-(^)%in2^^~2^1-(^)' ^ ^ 

where we restored the speed of light c and Q denotes the angular velocity of the pulsar. 
Working with this expression for the pulsar charge density, the resulting charge-separation 
can be determined as depicted in FIG. 3. By contrast, our force-free (and fully general rel- 
ativistic) model may be referred to as the solution of quasi-neutral plasma in which two 
species of charges can coexist at a given point in space while allowing for still non-zero net 
space charge density. And non-trivial poloidal current can still exist as it can be defined in 
terms of the difference in velocities between the two charge species. Then the origin/source 
of charges in this force-free treatment is mainly the pair creation due to strong electromag- 
netic field in space which provides infinite supply of ample plasma and hence the continuous 
flow of poloidal current without end. This difference in the nature of charges and the re- 
sulting poloidal currents between the two model is indeed the key to understand why our 
force-free and general relativistic treatment of the pulsar magnetosphere presents more self- 
consistent and hence satisfying view of the pulsar spin-down mechanism in the sense that 
it is consistent with the mechanism of Blandford and Znajek ^| employing the rotating 
black hole magnetosphere. This essential difference between the two pulsar models, how- 
ever, do not necessarily mean that our force- free and general relativistic treatment is able 
to provide a successful mechanism for pulsar spin-down in terms of the magnetic braking 
torque while Goldreich- Julian's original but non-general relativistic (GR) one fails to do 
so. And obviously in order to be convinced that both of the two models can successfully 
provide the mechanism for the pulsar's magnetic spin-down, one needs to demonstrate that 
the directions of poloidal currents particularly at the surface of the neutron star are indeed 
the same correctly leading to the magnetic braking torque that we discussed earlier. Thus 
in the following, this last point shall be addressed. Indeed, the directions of the poloidal 
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current (that closes globally in the magnetosphere) are defined differently in the two models. 
First in our force-free treatment, the direction of the poloidal current and the structure of 
the charge densities given in eqs.(21) and (29) respectively are determined simultaneously 
via the behavior of the quantity {dl/d^) as given in eqs.(35) or (36) just as it was the 
case with the rotating black hole magnetosphere 0, • Namely, one is not determined as 
a result of the other. Of course, this is because the continuous supply of ample plasma is 
assumed in our treatment as discussed above. And as depicted in FIG.l, the direction of 
poloidal current such determined particularly at the surface of the star correctly leads to 
the magnetic braking torque directed opposite to the star's angular momentum. 
In the original model of Goldreich- Julian's , on the other hand, the structure (i.e., 

the sign) of charge density is actually determined as a result of the direction of the poloidal 
current. That is, the strong local electric field near the surface of the neutron star first de- 
termines the direction of local poloidal current (i.e., the fiow of local charge carriers) which, 
in turn, determines the sign of local charges in a given domain of space. To be more specific, 
Goldreich and Julian began their analysis by assuming that the neutron star with an aligned 
dipole magnetic field is surrounded first by vacuum. Then the longitudinal electric field at 
the neutron star surface turns out to have component parallel to the poloidal magnetic field 
given hy E\\ ~ — cos^^ and particularly at the equator the vacuum electric field is directed 
radially outward. Thus in the polar region, the vacuum longitudinal electric field drives 
a current towards the star (by pulling space ions if present or by ripping electrons off the 
pulsar surface) while at the equator it drives current away from the star (by pulling space 
electrons if present or by ripping ions off the surface). In this way, the vacuum longitudinal 
electric field causes charge emission, i.e., the fiow of the poloidal current, only until the 
magnetosphere is filled with plasma with the charge density being given by the Goldreich- 
Julian's equilibrium value given in eq.(38). Namely, the poloidal current cannot fiow without 
end because once the charges get accumulated up to GJ's equilibrium value in eq.(38), the 
emission of charges becomes electrostatically unfavorable. And if, particularly in the particle 
acceleration region, there appears the difference between the local plasma charge density and 
the equilibrium Goldreich- Julian density p'^^ given in eq.(38), the longitudinal electric field 
arises and as a result, the plasma in the magnetosphere would be streaming out along open 
magnetic field lines past the light cylinder as a centrifugally slung, relativistic wind leading 
eventually to the observed radio emission. Note that if the current driven by the vacuum 
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longitudinal electric field can close in a global current system, particularly the direction of 
the current flowing on the stellar surface again would exert the correct magnetic braking or 
spin-down torque on the neutron star. In this non-GR Goldreich- Julian model, therefore, 
the charge-separation depicted in FIG. 3 does not really conflict with the spin-down process 
as the direction of the poloidal current is indeed the same as that of our force-free treatment 
discussed above correctly leading to essentially the same magnetic braking torque. Thus to 
summarize, regardless of the difference in the origin/source of the space charges and in the 
definition of the direction of poloidal current, both the Goldreich- Julian's original model 
and our present force-free treatment provide the working mechanism for magnetic pulsar 
spin-down. Nevertheless, we would like to stress again that there indeed is a more desirable 
feature in our treatment that distinguishes it from the original model of Goldreich and Ju- 
lian's. And it is the fact that the force-free and fully general relativistic treatment of the 
problem of pulsar magnetosphere presented in this work appears to provide much upgraded 
and closer view of the pulsar spin-down mechanism in the sense that it is consistent with 

12, Q| employing the rotating black hole 



the mechanism of Blandford and Znajek 
magnetosphere. 

Thus far, we have been interested in the comparison between our force-free and general rel- 
ativistic model and the Goldreich- Julian's non-GR model of pulsar magnetosphere. And it 
has been realized that the main differences between the two arise not from the GR effect but 
from the different nature and source of the charges. Now this realization leads us to turn to 
another relevant comparison. That is, it seems equally relevant to consider the comparison 
of our force-free and general relativistic treatment with a force-free but non-GR treatment 
of pulsar magnetosphere and see if there is actually something generic in fully general rela- 
tivistic treatment of the pulsar electrodynamics. Indeed such force-free but non-GR study of 
pulsar magnetosphere has been perfomed some time ago by Okamoto Sj and by Contopou- 



los, Kazanas and Fendt |23i]. Later on in the subsection 5 of the present section, the rigorous 
comparison of our present treatment with this second class of study shall be carried out and 
if we mention the essential result in advance, as far as the pulsar electrodynamics goes, the 
GR treatment does not seem to have any generic effect other than the stereotypical compli- 
cations and elaborations such as the large redshift near the neutron star's surface and the 
frame-dragging effect and hence the quest for the introduction of ZAMO, the local inertial 
observer carrying out the actual observations. Indeed, the existence of and the observations 
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by ZAMO is a non-trivial deviation from non-GR treatment of the pulsar electrodynamics 
since it is the strong electric field felt by ZAMO that actually renders the pair creation 
of charges via the so-called Schwinger process work [24]. Recall that our model of pulsar 
magnetosphere depends, for the source of ample supply of quasi-neutral plasma, heavily on 
the pair production of charges in space. 

To summarize, it has been quite uneasy to accept that the two relativistic spinning compact 
objects of nearly the same species, the neutron star (i.e., pulsar) (based on the Goldreich- 
Julian model) and the black hole (based on, say, the Blandford-Znajek model) have gener- 
ically different structures of magnetospheres. And in our force-free and general relativistic 
pulsar model, we realized that it shares the same structure of singular surfaces of flows with 
that of original Goldreich- Julian model on the one hand and shares essentially the same 
structures of charge-separation and the poloidal current with those of rotating black hole 
0, Q on the other. 





FIG. 3: Pulsar magnetosphere - the Goldreich-Julian's aligned rotator model. 

4. The energy and the angular momentum flux 

In the above, we showed, in terms of the space charge-separation structure of the pulsar 
magnetosphere, that in the force-free case the longitudinal (or poloidal) currents circulating 
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in the neutron star magnetosphere leads to the magnetic braking torque that actually spins 
it down in a similar manner to the case with the Blandford-Znajek mechanism for the ex- 
traction of rotational energy from Kerr black holes. In this subsection, we shall demonstrate, 
in terms of the energy and the angular momentum flux at the surface of the neutron star, 
that this argument does indeed hold true. The general expression for the redshifted energy 
flux Se and the angular momentum flux about the axis of rotation Sl are given respectively 



by 



= -^[a(E X B) - c<;(E ■ m)E - uj{B ■ m)B + -^^(E^ + B^)m], 
All 2 

Sl = — [-(E-m)E- (B-m)B + i(E2 + B2)m]. 
Att 2 



(39) 



Since the toroidal component of the fluxes are irrelevant, we only need to consider the 
poloidal components 
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where n denotes the unit vector outer normal to the neutron star surface. Now note that 
when the spin J of the rotating neutron star and the magnetic field B are parallel, B± > 0, 
J > whereas when J and B are antiparallel, B± < 0, J < due to their definitions eqs.(16) 
and (19). Namely, the magnetic flux (and B) is defined to be positive /negative when it 
directs upward/ downward while the poloidal current is defined to be positive /negative when 
it directs downward/upward as we noted earlier. Thus one always has IB±^ > 0, and hence 
from eq.(41) above, we always have 
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Since the angular momentum and the energy flux going into the neutron star surface are all 
negative, this means that the rotating neutron star (i.e., the pulsar) experiences magnetic 
braking torque, namely spins-down and as a result, always loses part of its rotational energy 
(at the surface). 

5. Limit of vanishing general relativistic effects 

In earlier subsections, we have studied the detailed comparison between our force-free and 
general relativistic model and the Goldreich- Julian's non-GR model of pulsar magnetosphere. 
It then has been argued that the main differences between the two arise not from the GR 
effect but from the different nature and source of the charges. This could be checked in a 
rigorous manner if we erase the GR content in our force-free treatment of the pulsar mag- 
netosphere and see if these differences still remain. We also have turned to another equally 
relevant comparison. Namely, we have considered the comparison of our force-free and gen- 
eral relativistic treatment with a force-free but non-GR treatment of pulsar magnetosphere 
to see whether there is actually something generic in fully general relativistic treatment of 
the pulsar electrodynamics. Again, such a comparison would be made explicit if the GR 

brce-free but non-GR study of 



22, M and hence the result of 



component in our treatment is washed out. Besides, such 
pulsar magnetosphere has been perfomed in the literature 
the comparison can be directly tested. Therefore in this subsection, we shall reconsider our 
force-free and general relativistic model and take its particular limit of vanishing GR content 
for these purposes. 

Evidently, taking the limit of vanishing GR content would amount to replacing the curved 
Hartle-Thorne spacetime exterior to the rotating neutron star with the fiat spacetime while 
maintaining the force-free nature of pulsar electrodynamics. And technically, this is equiva- 
lent to setting all the parameters associated with the non-trivial curved spacetime structure, 
i.e., the mass M, angular momentum J and the mass quadrupole moment Q in the Hartle- 
Thorne metric for the neutron star, to zero in all the equations of the pulsar electrodynamics 
presented in the subsections 1 and 2 above. 

In the limit of vanishing GR content, apparently the exterior spacetime is the fiat Minkowski 
one and in the following we shall take the cylindrical coordinates {t, R, 0, z) in which the 
Minkowski metric is given by 

ds^ = -df + dR^ + ^ ^^2_ ^43^ 
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Then the role played by proper distance from the axis of (neutron star's) rotation w that 
has been employed thus far in the general relativistic treatment shall henceforth be taken 
over by the radial coordinate R. Next, we start with the Maxwell equations in this flat 

spacetime 

V • E = 47rpe, V • B = 0, 
V X E = 0, V X B = 47rj (44) 

where we dropped the terms d{...)/dt — 0, d{...)/d(f) — due to stationarity and axisymme- 
try. Next, throughout, the force- free condition is still assumed to hold, i.e., 

PeE+jxB = 0, B = Bt + Bp (45) 

which also implies E ■ B = 0. Then this force- free condition indicates that the charged 
particles are sliding along the magnetic field lines and hence the toroidal (angular) velocity 
of magnetic field lines (which are frozen into plasma) is given by 

vf = R^pe^ = ftpm (46) 

and then j = Jt + jp with Jt = Pe^T where \t consists of given above and the streaming 
velocity along the toroidal magnetic field lines. Here, again m = Re^ = {OtptpY^'^e^- Then 
from the force-free condition above, it follows that 

E = Ep = -— Jt X Bp. (47) 

Pe 

We have established the force-free condition and based on this, we can now derive all the 
force-free pulsar electro dymanics equations. First we consider the poloidal field components. 

As before, a magnetic fiux through an area A whose boundary is a m-loop is given by 



= / B-dS. (48) 

J A 



Then from = • dr and alternatively ci^ = B • (dr x 2TrRe^) = {2'KRe^ x B) • dr, we 
get 

^''-^h^-^h^ ^^^^ 
where we used m • m = g^^^ = R^ and hence 

Ep = -vp X Bp = -^V*. (50) 
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We are now ready to discuss the poloidal current and the associated toroidal magnetic field. 
Once again, a poloidal current through the same area whose boundary is a m-loop is given 

by 



- [ j-dS (51) 

J A 



then similarly to the case of poloidal magnetic field, V/ = — 27rit!e ? x jp, we get 



. _ Vi><e^_ V/xm 

^""-'^^-'^^ ^^^^ 



and then from eqs(49), (52), 



jp = -^Bp. (53) 



Next, we turn to the Ampere's law in Maxwell equations 

J^{V X B) ■ dS ^ An Jj ■ dS (54) 

which, upon using the Stoke's theorem, becomes 27ri?|BT| = —Anl and hence 

X. 2/ 2/ 

Bt = -^e^ = -^m. (55) 
Next, using eqs.(44) and (50), one gets the Gauss law equation 

while the Ampere's law in eq.(44) yields 

,. = jL(VxB),^-Av.(i,V*). (57) 

Alternatively, by combining eqs.(56) and (57), one gets 

Pe = mpjr - ^B,. (58) 

We are now in a position to write down the force-free limit of the Grad-Shafranov equation. 
Take the force-free condition p^E -|- j x B = 0, and focus on its "poloidal component" 
equation, 

-Pe{VF X B)|p+j X B|p = 
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which gives, 



p,mF+3T + ^BT = Q. (59) 



Now by plugging eqs.(56) and (57) in (59) above and using eq.(55), one arrives at 



V- 



IQn^I dl , , 



This is the stream equation that, in principle, would allow us to determine the field struc- 
ture of the force- free pulsar magnetosphere. Lastly, from eqs.(58) and (59), one gets the 
expressions for the charge and toroidal current density 

where we restored the speed of light c for the sake of comparison with their counterparts in 
the Goldreich- Julian model and as usual Rl = c/Qp denotes the radius of the light cylinder. 
We now test our force-free and general relativistic model by comparing its limit of vanishing 
GR content we have studied in this subsection firstly with the force-free but non-GR model 
of [2^ [2^ in (I) and then next with non-GR model of Goldreich and Julian's in (II) below. 
(I) Among other things, it is noteworthy that all the pulsar electrodynamics equations and 
particularly these expressions in eq.(61), except for the simplifications due to the absence 
of the GR content, remain essentially the same as their fully GR counterparts given earlier 
in the subsections 1 and 2. This implies that the GR treatment does not seem to have 
any generic effect on the pulsar electrodynamics other than the stereotypical complications 
and elaborations such as the large redshift near the neutron star's surface and the frame- 
dragging effect and hence the quest for the introduction of ZAMO, the local inertial observer 
carrying out the actual observations. Thus in a sense, our present model can be thought of 
as a formal general relativis tic gene ralization of the force-free but non-GR model of pulsar 



magnetosphere suggested in 



23|. Notice that the expressions in eq.(61) above essentially 
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23|. 



coincide with the corresponding results constructed in 
(II) Next, one can immediately realize that the equilibrium GJ charge density given in 
eq.(38) is just a special (vacuum space charge) case of eq.(61) in which / = 0. Note also 
that the charge density in eq.(61) above in our force- free treatment actually can be written 



as 



Q 4ttI dl 

pr = ^ ^V\2 +p^' (62) 
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with p^-^ being the equihbrium GJ charge density given in eq.(38). Recall here that pf"^ 
represents the maximum amount of charge available in the Goldreich- Julian's vacuum pulsar 
model and the origin/source of these charges is basically the surface of the star itself. In 
our force- free model, however, the origin/source of charges is mainly the pair creation due 
to strong field in space and hence it presumably guarantees the infinite supply of ample 
plasma and hence the continuous ample flow of poloidal current. Next, determining the 
structure of charge separation (i.e., the sign of separated charges) in our force- free pulsar 
model using eq.(61) does not look so simple as we need to figure out which term, between 
\{ATT/c^)I{dI /d^)\ and \Bz\ in the numerator of eq.(61), is greater than the other in different 
domains of the magnetosphcrc. As wc mentioned earlier, fortunately there is a guiding 
principle that allows us to determine the sign of separated charges in different domains 
with confidence and it is the insightful realization that with respect to ZAMO, the local 
inertial observer, the magnetosphere around the rotating neutron star should rotate in the 
same direction as the compact object itself. Namely, using the definition of the angular 
momentum of an electromagnetic field, Jem = / d^x[r x (Ep x Bp)] and demanding Jem ~ 
Jns, one can determine in an unambiguous manner the structure of the charge separation 
in different domains of the pulsar magnetosphere as summarized in FIG.l and it turns out 
to be essentially the same as that in the GJ model (but only inside the light cylinder) given 
earlier in FIG. 3. 

To summarize, it should now be clear that the main differences between the two models 
(i.e., ours versus GJ's) arise not from the GR effect but from the different nature (such as 
the force-free assumption) and source of the charges. But the essential features such as the 
structure of charge separation and the direction of the poloidal current (particularly at the 
pulsar surface) leading to the pulsar spin-down due to the magnetic breaking are shared by 
the two models. Next, it seems worth contrasting carefully the nature of the critical field 
lines in the two pulsar models. First, the critical magnetic surface in our model is the surface 
on which (dl/d^) — 0. On the other hand, in GJ model of purely charge-separated pulsar 
magnetosphere, the critical field line has been defined as the one on which flpBz = (see 
eq.(38)). Thus the critical magnetic surface ^ = \l'c in these two models generally may not 
coincide. Indeed, on the critical magnetic surface \E' = in our force- free treatment, the 
net current is zero (i.e., inflowing current = outflowing current) but the charge density there 
may not be exactly zero and hence the critical magnetic surface in FIG.l may not serve 
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as exactly a "charge-separating" surface. In the simpler model of Goldreich and Julian's, 
however, the critical field line in FIG. 3 is indeed precisely a charge-separating boundary. 

IV. Pulsetr jet equation - The general Grad-Shafranov equation 

In this more general Grad-Shafranov equation, the role played by the plasma particles, 
i.e., their dynamics, has been taken into account. 
1. Basic equations 

First, the Maxwell equations in the background of the stationary axisymmetric rotating 
neutron star spacetime given earlier in eq.(12) should be supplemented by the charge con- 
servation 

V^J^ = 0, or ^ + V-(aje) = 0. (63) 

The remaining general relativistic magnetohydrodynamics (MHD) equations are ; 

{Particle (mass) conservation) 

d 

Va(nw") = 0, or ^(7^) + V • (a-fnv) = (64) 

where = dx^/dr — {cj,jv) denotes the fluid 4- velocity and UaU" — gapu^u^^ — — 1, 
7 = (l-vVc2)-V2. 

{Energy-momentum conservation) 



Tf = (^) + Pq'-'^ (65) 
= ^[F;F^^--g'^f{F,^F^^^)] such that 



47r 

where w — {e-\- P)/n is, the specific enthalpy in which P denotes the proper pressure and 
e denotes the proper internal energy density given by e = nmc^ -|- (F — 1)~^P and hence 
w^mc^ + VP[n{V -l)]-^. 
{Infinite conductivity (Ideal MHD)) 

F^uf^ = 0, or E + -V X B = 0. (66) 
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{Equation of state (Entropy conservation)) 



s{w,P)^kB{r-l)-Hn{Pn-^) (67) 

where F = 5/3, 4/3 for non-relativistic motion and for ultrarelativistic motion, respectively. 
Then by contracting Ua with eqs.(65) and (66) and using the 1st law of thermodynamics 
dw — Tds + -dP, one gets 

d 

Va{nsu") — Q, or — (7ns) + V • (Q;7nsv) = 0. (68) 

{Momentum conservation (Euler equation)) 

By contracting the energy-momentum conservation equation (65) with (gf^A + UaU\) and 
then employing the Maxwell equations, one gets 

nw{uNp)ua = -daP - Ua{u^Vp)P + ^F^^ Jf . (69) 

Particularly in the "cold hmit"(P = 0, e = nmc^, and w — e/n — mc^), it reduces to 

nmc^uf^V p)u^ = -F^pJ^. (70) 



2. The Grad-Shafranov (GS) approach 

In this section, we are mainly interested in the derivation of the Grad-Shafranov (GS) 
equation which describes the dynamics of plasma particles. And in the following, all the 
time derivative terms will be dropped, i.e., d{...)/dt = due to the stationarity of the 
background Hartle-Thorne metric for the region exterior to the rotating neutron star. 
2.1 Constants of motion 

(I) Substituting E = Ep = - 
one can readily realize that 

B • VQf = (71) 

indicating that ilp is constant on magnetic surfaces, i.e., flp — Of(^) which represents the 
generahzed Ferraro's isorotation. 

(II) Combining 

the freezing-in condition ; E^ -|- ^(v x B)t = 0, 

32 



2na 



V* into the Maxwell eq.(12) V x (aE) = (B • Vu;)m, 



the particle conservation ; V • (a'jnv) — 0, 

and the Maxwell equation ; V • B = 

one ends up with Up = = q (Bp/an) and hence from 



Ur = 7Vt = r] f — By ) + 7 
Van 



Vtp — uj' 



a 



wet 



it follows that 



u — ^v — — B + 7 

an 



rip — LV 



a 



we- 



(72) 



(73) 



where the quantity represents the particle flow along the magnetic flux or the particle-to- 
magnetic field flux ratio. 

Then plugging (73) back into the particle number consevation eq.(64) yields 



= V • (anu) = V ■ (r/B) 

= 77(V • B) + B • {Vf]) = B • (V77) 



(74) 



which implies that 77 must be constant on magnetic surfaces as well, i.e., rj — r){^). 

mm 

Let x'^ be a Killing field associated with an isometry of the background spacetime metric, 
then 



which yields V^(T^'^Xm) = 0- 



(75) 



Since the Hartle-Thorne metric possesses the time-translational isometry and the rotational 
isometry, there are corresponding Killing fields /c'* = (d/dt)^ and m'^ — {d/d(j))^ respectively, 
such that the quantities 



ei'^-Ti'^K and ^T^^'m^ 
are covariantly conserved. To be a little more precise. 

Thus using, 

( V 

— (c7, 7v) = C7. — B + 7 

\ na 



(76) 



(77) 







a 





(78) 
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and 



-To" = nu'E, 



(79) 



C = t; = nu'L 
one gets two more integrals of motion P| 



E = E{^) 
L 



2tt 
I 



(80) 



L(^) = — + wrfwu^ 



and the total loss of energy and angular momentum are given by 

m n.T. 



tot 







tot 



^ max 



E(^)d^, 
L(*)(i*. 



(81) 



(V) The entropy conservation Va(rasu°) = reduces, for stationary axisymmetric case, to 



V ■ iansyy] = 0. 



Thus using 



u 



an 



f2p — U! 



a 



roe?, 



one gets 



= V ■ (ansu) = V ■ {r]sB) 

= sV ■ (r/B) + r]B ■ (Vs) = r]B ■ (Vs) 



(82) 



(83) 



(84) 



which implies that the entropy per particle s must be constant on magnetic surfaces as well 



s = si"^). 



(85) 



To summarize, for the stationary axisymmetric case, there are 5-integrals of motion (con- 
stants on magnetic surfaces) 



{^^(vl/), r/(v[/), s(vl>), E(v[/), L(vl/)}. 



(86) 



We shall now show that once the poloidal magnetic field Bp and the 5-integrals of motion 
given above are known, the toroidal magnetic field i?<i and all the other plasma parameters 
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characterizing a plasma flow can be determined. 

To do so, we solve the two conservation laws in eq.(80) and the toroidal component of eq.(73) 



V 

= — + 7 



for {/, 7, u^} to get 



an 



fip — Co" 



a 



zu 



2r]I 
a'^nzu 



7 



a 



zu 



(87) 



I _ a'^L-{VLF-uj)w'^{E-ujL) 
2^ ~ 
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1 a^E-npL)- M^{E-ujL) 
arjw a^ — iyip — loY"^'^ — M"^ ' 
1 {E - VLpL) {VLf - iu)zu^ - LM^ 



(88) 



zur]w a^ — {^p — Luyzu"^ — 

where = ATirfw/n = a'^{Up/u\) is the square of the Mach number of the poloidal velocity 
Up = r](Bp/na) with respect to the Alfven velocity ua = Bp^Annw)'^^"^. 
Now in order to determine this Mach number, consider 



2 2 2 2 2 2lA 2\ 

7 — u =7 —7V =7(1 — v) = 

and into this relation, we substitute eqs.(73) and (88) to get 

K 1 M^(V*)2 



(89) 



t^M2 647r4 



I 2 2 2 

^ a 1] w 



(90) 



where 



A = a^ -{Qf- ojf^"^ -M^ = N'^ - M^, 
K = a^w'^{E - VLpLfia^ - {VLp - ujfw'^ - 2M'^] 
+ M^[w\E - uLf - a^L^] 

which is the Bernoulli equation. 

To summarize, once Bp, {Qpi"^), vi'^)^ -^l^)) -^(^)) -^(^)} known, the characteristics 
of the plasma flow, {/(or i?^), 7, u^, Up, M^(or ua)} can be determined by eqs.(88)-(90). 
2.2 The Grad-Shafranov equation 

The Grad-Shafranov equation is the "trans-fleld" equation of magnetic fleld lines and it 
results from the poloidal component of the Euler equation (69). Further the Grad-Shafranov 
equation describes a "force-balance" in the transfleld (i.e., poloidal) directions. For the case 
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at hand in which the content of jjlasma dynamics is taken into account, the Grad-Shafranov 
or the pulsar jet equation reads 



n 1 dn o ^ ds 

r]d^ d\I> 

where T denotes the temperature and G = a^w'^{E - VLpLf + a^M'^L'^ - M'^w'^{E - uoLf. 
Note that this Grad-Shafranov equation contains only \E' and 5-integrals of motion, position 
and physical constants. Thus the Grad-Shafranov equation is autonomous. 
Also it is interesting to note that taking the limit, and s — > 0, this pulsar jet 

equation given above reduces to the pulsar equation, i.e., the force-free limit of the Grad- 
Shafranov equation (neglecting the content of plasma dynamics) given earlier in eq.(28) 
and at the same time the 5-integrals of motion also reduce to just 2-integrals of motion 
(ili?(\E'), which can be envisaged from eq.(80). 

2.3 Singular surfaces 

The algebraic equations (88) and (90) allow for the determination of the locations of the 
singular surfaces of general relativistic MHD flows. 

{Alfven surfaces) 

From eqs.(88) and (90), one realizes that there exists general relativistic version of the Alfven 
points where A = a'^ — {Q,p — ijjfw'^ — = holds. Then using = ^^(u^/m^)^ one 

3 



immediately sees that on the Alfven surface 



,2 „,2 



2 h (^F -^)^^^ ' 

UZ = Ua 11 



(92) 



must hold which, in the non- relativistic limit, coincides with the Alfven velocity. On this 
Alfven surface, in order to keep the value of {/, 7, m,^} in eq.(88) finite, one requires that 
numerators vanish there as well. This constraint amounts to a single relation [l7| 

+ uj{ilF - oj)vj^ L - {VLf - uj)w'^E = 
or equivalently ilp [L E) = — — — -. (93) 

Note that it possesses essentially the same structure as its (rotating) black hole counterpart 



17[ . This is a general relativistic version of the Newtonian result that the angular momentum 
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carried away by the wind is given by the position of the Alfven point [25]. Eqs.(92) and (93) 
also allows us to express the location of a single Alfven point as 

1/2 



ZUa 



(94) 



{np - uj){E - ioL) 
{Light cylinders) 

Like in the force-free case we discussed earlier, the pulsar magnetosphere under consideration 
possesses a single light cylinder whose location is given by A^^ = — [Qp — uYw'^ = 0, 
namely at 

as fii? > everywhere for the case of rotating neutron star as we stressed earlier. And in the 
force-free limit, and s — or equivalently E = QpL, the Alfven surface discussed 

above coincides with this light cylinder, i.e., zua = ^l- Next, the possible existence of the 
fast and the slow magnetosonic surfaces in this case of pulsar magnetosphere can be checked 
following essentially the same procedure as that in the case of rotating (Kerr) black hole 
magnetosphere. Perhaps, the easiest way of defining these magnetosonic surfaces is to think 
of them as being singularities in the expression for the gradient of the Mach number M. 
Here, however, we shall not go into any more detail and instead, we refer the interested 
reader to B and Q for related discussions. 
{Injection surfaces) 

Lastly, we introduce the injection surfaces, r = rj[6,QF{'$)] for both plasma inflow and 
outflow where a poloidal flow starts with a sub-Alfvenic velocity. And the plasma inflow or 
outflow which starts from this injection point must pass through the Alfvenic point to reach 
the neutron star surface or the far region. In order to determine these surfaces, however, we 
need some concrete physical model which is beyond the scope of the present work. 
3. Problems with the Grad-Shafranov approach 

We now discuss the difficulties when treating the (rotating) black hole or pulsar magneto- 
sphere in terms of the so-called Grad-Shafranov approach. As has been pointed out thus 
far, the central role is played by the Grad-Shafranov equation in determining the structure 
of electromagnetic field and the characteristics of the plasma flow in the black hole or pulsar 
magnetosphere. Thus we begin with the algorithm to solve the Grad-Shafranov equation, 
(i) Once the physical constants {n, w, T, Bp} are known and the 5-integrals of motion 
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{rii;'(^f), ?7(^f), s(^), E{'^), are given, 

(ii) one might be able to solve the Grad-Shafranov equation in eq.(91) for the poloidal mag- 
netic flux or the stream function \1/ = \l/(r, 6) as a function of the poloidal coordinates (r, 6). 

(iii) Then from this ^ = ^(r, 6) and using 

tip = — , tiT 



E 



271ZU azu 



p 



V^, Et = (96) 



27ra; 

one in principle determines the structure of the electromagnetic fields and then next using 
eqs.(88)-(90), one obtains the characteristics of the plasma flow {/, 7, u^, Up, M^}. 
In this way, in principle, one can determine the structure of pulsar/black hole magnetosphere. 
In practice, however, this Grad-Shafranov approach does not appear to be so tractable since 
in the step (i), there is no known systematic way of evaluating the "physical constants" and 
giving the "5-integrals of motion" in terms of the stream function 

In the force-free case we discussed earlier, however, the plasma content is now absent and the 
whole task of dealing with the Grad-Shafranov approach reduces to the attempt at finding 
the solution (i.e., the stream function \l/(r, ^)) of the stream equation (28). Even in this 
simpler case, one is still left with the ambiguity in determining the 2-integrals of motion, 
{fii;'(\l/), /(^I/)} in a self-consistent manner. Indeed, it is instructive to note that the stream 
equation (28) is nonlinear but the nonlinearity entirely comes from the integrals of motion. 
Thus in the simplest, non-realistic case when the current /(^P) is absent and the field angular 
velo city is constant fii?(\I') = Q, the stream equation (28) becomes linear and thus soluble 
|26l l27|. Then one might wish to elaborate on this simplest case to construct more general. 



realistic solutions by "guessing" consistent current ansatz I = /(\1/). Such attem pts act ually 
have been made and for more details in this direction, we refer the reader to 

V. Summary and discussion 



em pts actually 



In the present work, we performed the study of the pulsar magnetosphere in the con- 
text of general relativistic magnetohydrodynamics (MHD) by employing the so-called Grad- 
Shafranov approach. We considered both the force-free and full MHD situations and ac- 
cordingly derived the pulsar equation and the pulsar jet equation respectively. The resulting 
Grad-Shafranov equations and all the other related force-free equations or general relativistic 
MHD equations turn out to take essentially the same structures as those for the (rotating) 

38 



black hole magnetosphere. The essential distinction between the two cases, however, is the 
spacetime (metric) contents. For the pulsar magnetosphere case, one needs to choose the 
Hartle-Thorne metric mentioned above whereas for the black hole magnetosphere case, one 
has to select the Kerr black hole metric. In this way, we demonstrated that the pulsar and 
the black hole magnetospheres can be described in an unified and consistent manner. 
Also there is quite an uncomfortable state of affair that there has been no complete model 
for the structure of longitudinal (or poloidal) currents circulating in the neutron star mag- 
netosphere that can provide the solution to the problem, say, of pulsar spin-down. To this 
problem, we have provided a partly satisfying solution again by treating the region outside 
a magnetized rotating neutron star as a curved spacetime represented by the Hartle-Thorne 
metric. Namely, we have demonstrated that both for pulsar and for rotating black hole 
cases the structure of charge-separation and the direction of longitudinal (or poloidal) cur- 
rent circulating in the magnetospheres actually lead to the magnetic braking torques that 
spin down the rotating neutron star and the black hole regardless of whether the spin and the 
(asymptotic) direction of the magnetic field are parallel or antiparallel. We also remarked 
that the structure of charge-separation that resulted from our force-free treatment of the 
pulsar magnetosphere turns out to be the same as that in the original model of Goldreich 
and Julian jll]. And this unified picture can be thought of as a more satisfying solution to 
both the magnetized rotating neutron star interpretation of radio/X-ray pulsars and the 
rotating supermassive black hole interpretation of AGNs/quasars or even GRBs [l^ . 
Next, one might be worried about the validity of the Hartle-Thorne metric for the region 
surrounding the slowly-rotating neutron stars employed in this work to describe the mag- 
netosphere of pulsars which seem rapidly-rotating having typically millisecond pulsation 
periods. Thus in the following, we shall defend this point in a careful manner. Here the 
"slowly-rotating" means that the neutron star rotates relatively slowly compared to the 
equal mass Kerr black hole which can rotate arbitrarily rapidly up to the maximal rotation 
J = M^. Thus this does not necessarily mean that the Hartle-Thorne metric for slowly- 
rotating neutron stars cannot properly describe the millisecond pulsars. To see this, note 
that according to the Hartle-Thorne metric, the angular speed of a rotating neutron star 
is given by the Lense-Thirring precession angular velocity in eq.(9) at the surface of the 
neutron star, which, restoring the fundamental constants to get back to the gaussian unit. 
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IS 




with J = aM 



If) (°<^<^) 



(97) 



As we mentioned earlier, one of the obvious differences between the black hole case and the 
neutron star case is the fact that the black hole is characterized by its event horizon while 
the neutron star has a hard surface. As such, in terms of the spacetime metric generated by 
each of them, just as the Lense-Thirring precession angular velocity (due to frame-dragging) 
at the horizon represents the black hole angular velocity, the Lense-Thirring precession an- 
gular velocity at the location of neutron star's surface should give the angular velocity of 
the rotating neutron star. 

Thus the Hartle-Thorne metric gives the angular speed of a rotating neutron star, hav- 
ing the data of a typical radio pulsar, M ~ 2 x lO^^(g'), ro ~ lO^(cm), as a; = 
2a{M'^ /rl){G/c){G/c^) ~ 10^(l/sec) which, in turn, yields the rotation period of r = 
27r/a; ~ lQ-^{sec). And here we used, (G/c^) = 0.7425 x lQ-'^^{cm/g) and (G/c) = 
2.226 X l{)~^^{cw? /g ■ sec). Indeed, this is impressively comparable to the observed pul- 
sation periods of radio pulsars r ~ 10~^ — l(sec) we discussed earlier. As a result, we expect 
that the Hartle-Thorne metric is well-qualified to describe the geometries of millisecond pul- 
sars. 

Lastly, although the Grad-Shafranov approach toward the study of the pulsar magneto- 
sphere is not fully satisfying for reasons stated earlier, it nevertheless is our hope that at 
least here we have taken one step closer toward the systematic general relativistic study 
of the electrodynamics in the region close to the rotating neutron stars in association with 
their pulsar interpretation. 
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